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ABSTRACT 

In the present work, motivated by the work of Cai and Su [Phys. Rev. D 81, 103514 (2010)], 
we propose a new type of interaction in dark sector, which can change its sign when our universe 
changes from deceleration to acceleration. We consider the cosmological evolution of quintessence 
r*j ■ and phantom with this type of interaction, and find that there are some scaling attractors which can 

help to alleviate the cosmological coincidence problem. Our results also show that this new type of 
interaction can bring new features to cosmology. 
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I. INTRODUCTION 

In dark energy cosmology (see e.g. [l[ for reviews), the cosmological coincidence problem is one of the 
well-known conundrums, which asks: why are we living in an epoch in which the densities of dark energy 
and matter are comparable? To alleviate this coincidence problem, it is natural to consider the possible 
interaction between dark energy and dark matter in the literature (see e.g. pj-Ull ) . In fact, since the 
nature of both dark energy and dark matter is still unknown, there is no physical argument to exclude 
the possible interaction between them. On the contrary, some observational evidences of the interaction 
in dark sector have been found recently. For instance, Bertolami et at [12| showed that the Abell Cluster 
A586 exhibits evidence of the interaction between dark energy and dark matter, and they argued that this 
interaction might imply a violation of the equivalence principle. On the other hand, Abdalla et al. [l3j 
found the signature of interaction between dark energy and dark matter by using optical, X-ray and weak 
lensing data from the relaxed galaxy clusters. So, it is reasonable to consider the interaction between 
dark energy and dark matter in cosmology. 

In the literature, it is usual to assume that dark energy and dark matter interact through a coupling 
term Q, according to 

p m + 3Hp m = Q , (1) 
p de + 3H(p de + p de ) = -Q , (2) 

where p m and pd e are the densities of dark matter and dark energy (we assume that the baryon component 
can be ignored); pde is the pressure of dark energy; H = a /a is the Hubble parameter; a is the scale 
factor; a dot denotes the derivative with respect to cosmic time t. Notice that Eqs. ((TJ and © preserve 
the total energy conservation equation p t ot + 3H(p to t + Ptot) — 0, where ptot — Pm + Pde- Since there 
is no natural guidance from fundamental physics on the interaction Q, one can only discuss it to a 
phenomenological level. The familiar interactions extensively considered in the literature (see e.g. 0-[Hl) 
include Q — 3aHp m , Q — 3j3Hp to t, and Q = 3i]Hpd e - 

Recently, Cai and Su [14| investigated the interaction in a way independent of specific interaction 
forms by using the latest observational data. They divided the whole range of redshift z into a few 
bins and set the interaction term 8(z) = Q/(3H) to be a constant in each redshift bin. From the latest 
observational data, they found that 5{z) is likely to cross the non- interacting line (5 = 0), namely, the 
sign of interaction Q changed in the approximate redshift range of 0.45 <z< 0.9. In fact, this result 
raises a remarkable problem. Indeed, most interactions extensively considered in the literature, such as 
Q = 3aHp m , Q = 3 (5 H ptot and Q = 3i]Hpd e , are always positive or negative and hence cannot give the 
possibility to change their signs. As noted by the authors of [l4j |. some new interaction forms should be 
proposed to address this problem. 

In the present work, we are interested to propose such a type of interaction and consider its implications 
to cosmology. The authors of [lij found that the sign of interaction Q changed in the approximate redshift 
range of 0.45 < z < 0.9. We note that this redshift range is coincident with the one of our universe changing 
from deceleration to acceleration So, a simple idea naturally comes to our mind. If the interaction Q 
is proportional to the deceleration parameter 

a H 

= -l-7H^-!> (3) 



a 



H 2 H 2 



the sign of Q can change when our universe changes from deceleration (q > 0) to acceleration (q < 0). 
Noting that the deceleration parameter q is dimensionless, from Eqs. ((T|) and ©, Q oc qp and Q oc qHp 
are both viable from the dimensional point of view. To be general, we consider the linear combination of 
these two, namely 

Q = q(ap + 3pHp) , (4) 

where a and (3 are both dimensionless constants. It is not surprising to find p in the interaction Q. 
We note that in the literature (see e.g. [l(|) the derivatives of energy density have already been allowed 
to appear in the general forms of Q. However, the key point of our new interaction is the deceleration 
parameter q in Q, which makes our proposal different from the previous works. This new feature gives the 
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possibility that interaction Q can change its sign, and hence brings some interesting results to cosmology. 
In this work, we would like to consider three interactions of this type, namely 

Q = q(ap m + 3/3Hp m ) , (5) 

Q = q{aptot + SPHptot) , (6) 
Q = q(ap de + 3/3Hp de ) . (7) 

In the present work, we consider the cosmological evolution of quintessence and phantom with the 
above type of interaction. In Sec. [TTl we present the dynamical system of interacting quintessence and 
phantom. In Sees. IIIII — TVl we discuss the cases with Q given in Eqs. JSJ) — 0, respectively. In Sec. IVI[ 
a brief conclusion is drawn. 



II. DYNAMICAL SYSTEM OF INTERACTING QUINTESSENCE AND PHANTOM 

In this work, we consider a flat Friedmann-Robertson- Walker (FRW) universe. The Friedmann and 
Raychaudhuri equations are given by 

K 2 K 2 

H 2 = Y Ptot = Y^ Pde + Pm ^ ' ^ 

2 2 
H = —{ptot + Plot) = 2~(An + Pde + Pde) , (9) 

where k 2 = 8ttG. The role of dark energy is played by quintessence or phantom, namely 

Pde =P^ \e4? + V(<t>) , (10) 

p d e=P<t, = \ej> 2 -V{<j>), (11) 

in which e = +1 (quintessence) or e = —1 (phantom); V(<p) is the potential. In this work, we consider 
the exponential potential 

V{(j>) = V e- XK ^ , (12) 

where A is a dimensionless constant. Without loss of generality, we choose A to be positive, since we can 
make it positive through field redefinition <\> — > —(f) if A is negative. 

We consider the cosmological evolution of interacting quintessence (phantom) by using the method of 
dynamical system [l5| . Following 043; j we introduce the following dimensionless variables 



k0 nyV kJ~p~^ 

y=—K77> z = TtTTt ■ ( 13 ) 



With the help of Eqs. ((8j — (|TT|) . the evolution equations ((T|) and ([2]) can then be rewritten as a dynamical 
system, namely 

x ' = (s-3)x+- e (^Xy 2 -Q i y (14) 

y' = sy - \j\ Xx v > ( 15 ) 

s - I) z + Q 2 , (16) 
where 

*Q_^ q (17) 
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a prime denotes derivative with respect to the so-called e-folding time N = In a, and 

s=--^=3ex 2 + ^z 2 . (18) 

The Friedmann constraint equation © becomes 

ex 2 + y 2 + z 2 = 1 . (19) 
The fractional energy densities of dark energy and dark matter are given by 

2 2 
n r> — K P<p 2,2 r>— K 2 /orA 

ft de = fi = = ex + y , Cl m = -j^- = z . (20) 

Once the interaction Q is given, we can obtain the critical points (5, y, z) of the autonomous system 
Eqs. (|14p — (TT5|) by imposing the conditions x' = y' = z' = 0. Of course, they are subject to the Friedmann 
constraint Eq. (|T5|) . namely, ex 2 + y 2 + z 2 = 1. Note that these critical points must satisfy the conditions 
y > and z > by definition (|13l) . and the requirement of 5, y, z all being real. Then, we can discuss 
the existence and stability of these critical points. An attractor is one of the stable critical points of the 
autonomous system. 

To study the stability of the critical points of Eqs. (fl"4|) — (fTB]) . we substitute the linear perturbations 
x —¥ x + dx, y — > y + 5y and z — ¥ z + Sz about the critical point (x, y, z) into Eqs. (|14[) — (|16l) and linearize 
them. Because of the Friedmann constraint Eq. (1191) . there are only two independent evolution equations, 
namely 

5x' = (s - 3)Sx + xSs + - (V6XySy - SQij , (21) 
Sy' = ySs + sSy - A (xSy + ySx) , (22) 

where 

s = (ex 2 - f + 1) , Ss = 3 {exSx - ySy) , (23) 

and SQi is the linear perturbation coming from Q\. The two eigenvalues of the coefficient matrix of the 
above equations determine the stability of the critical point. 

In the following sections, we will study the dynamics of quintessence (phantom) with the interaction 
Q given in Eqs. §5§ — (JTJ), respectively. 



III. THE CASE OF Q = q(ap m + 3/3Hp m ) 

Firstly, we consider the case of Q — q(ap m + 3f3Hp m ) given in Eq. (|5|). Substituting it into Eq. ([1]), 
one can find that 

Prn = 7^ ■ '!H Pm . (24) 
1 — aq 

Then, substituting into Eq. ((5]), we can finally obtain 



1 — aq 



3qH Pm . (25) 



At first glance, this interaction form is very similar to the familiar Q — 3i]Hp m in which 77 is a constant. 
However, the deceleration parameter q in Eq. (j25|) makes difference. Note that q = —l—H/H 2 is a variable 
function of time, which changes its sign when the universe changes from deceleration to acceleration. 
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Label 



Critical Point (x, y, z) 



M.lp 
M.lm 
M.2p 
M.2m 
M.3p 
M.3m 
M.4 
M.5 
M.6 



T2 



1)A: 



T2 



+ +r 2 - l)/e, 
+r 2 - l)/e, 



0, 
0, 

o, 

0. 

o, 

0. 



V2 - n + r 2 



V2 - n 
V2-ri 



ri + r 2 



- r 2 



- r 2 



A/(V6e). v/l-A 2 /(6e). 



/372(n -r 2 )/A, V 1 " n + ^ + 3 e (n - r 2 ) 2 /(2A 2 ) , - r 2 - 3 e (n - r 2 ) 2 /A 2 



/372(ri +r 2 )/A, ^/l - n - r 2 + 3 e (n + r 2 ) 2 /(2A 2 ) , ^/n + r 2 - 3 e (n + r 2 ) 2 /A 2 



TABLE I: Critical points for the case of Q = q(ap m + 3/3Hp m ) . 



Substituting Eq. (|25p into Eq. (jTT)) . we find that the corresponding Q% and Q 2 are given by 

3 (/3 - a)g z 2 3z (/3 - a)g , . 

Qi = o ' "1 ' = ^ : • (26 

2 1 — aq x 2 1 — a? 

Notice that q = s — 1 and s is given in Eq. (|18p . Then, substituting them into the autonomous system 
Eqs. (fl4|) — (|16p . we can find the critical points and present them in Table [I] Note that r\ and r2 are 
given by 



2 + 2a + 3/3 JAa 2 + (2 + 3/3) 2 - 4a(4 + 3/3) 

ri = ^ , »'2 = 7. • (27 

Qa oa 

If a = 0, only the first two critical points (M.lp) and (M.lm) in Table U can exist, which are trivial 
solutions in fact. If a ^ 0, for convenience, we can regard n and f2 as the model-parameters, in place of 
a and /3. By reversing Eq. (|2"T1) . we can express a and j3 as functions of ri and r 2 , namely 



2(2 - 6ri + 3r 2 - 3r ; 2 
4 - 12n + 9r 2 - 9r 2 ' ' 4 - 12n + 9r 2 - 9r 2 



-> ■ ^ /I _ 1 _L fw2 _ n^2 ■ ( 28 ) 



Now, we discuss the existence of the critical points in Table |TJ Obviously, Points (M.lp) and (M.lm) can 
exist only for e = +1 (namely quintessence). They are both quintessence-dominated solutions, because 
the corresponding fl m = z 1 = 0. For Points (M.2p) and (M.2m), if e = — 1 (namely phantom), we should 
have r\ — ri < 1 by requiring x to be real. However, in this case Q m — z 1 > 1 which is physically 
meaningless. Therefore, Points (M.2p) and (M.2m) can exist only for e = +1 (namely quintessence) and 
1 < r\ — r 2 < 2. Similarly, Points (M.3p) and (M.3m) can exist only for e = +1 (namely quintessence) 
and 1 < r% + r 2 < 2. For Point (M.4), if e = +1 (namely quintessence), it can exist under condition 
A 2 < 6; on the other hand, if e = — 1 (namely phantom), it can exist for any A. In fact, it is a dark- 
energy-dominated solution, because the corresponding il m — z 2 — 0. Point (M.5) exists under condition 
1 - n + r 2 + 3 e (n - r 2 ) 2 /(2A 2 ) > and 1 > r x - r 2 - 3 e (n - r 2 ) 2 /A 2 > 0. Point (M.6) exists under 
condition 1 - r 1 - r 2 + 3 e (n + r 2 ) 2 /(2A 2 ) > and 1 > n + r 2 - 3 e (n + r 2 ) 2 /A 2 > 0. Obviously, Points 
(M.2p), (M.2m), (M.3p), (M.3m), (M.5) and (M.6) are all scaling solutions, because the corresponding 



To study the stability of these critical points, by linearizing Q\, we obtain 

_ . AT 1 1 

(29) 



xn 3 P~ a J/l -2 -2\ 6 1 



2a; 1 — ag [ 1 — 



2exfa + 2ydy + (l - ex 2 - y 2 ) ^ 



where g = s — 1 and <5q = (5s, while s and 5s are given in Eq. (|23p . Substituting this <5Qi into Eqs. (|2ip 
and (|22D . the two eigenvalues of the coefficient matrix of Eqs. (l2"TT) and (|2^|) determine the stability of the 



G 



Point 


Eigenvalues 


M.lp 




-3(4- 


- 4n +ri- ri)/(-in + 3ri - 3rj{) , 3 - A^3/(2e) 




M.lm 




-3(4- 


- in + r? - rf)/(-4r a + 3r 2 - 3rf) , 3 + A x /3/(2e) 




M.2p 


6r 2 (- 


-2 + n -r 2 )/(2 


- 5n + 3r 2 + r 2 - 3rf) , (3/2) [n - r 2 - A^n - r 2 - 


l)/(3e)] 


M.2m 


6r 2 (- 


-2 + n -r 2 )/(2 


- 5n + 3r 2 + r 2 - 3rf) , (3/2) [n - r 2 + A v / 2(n - r 2 - 


l)/(3e)] 


M.3p 


6r 2 (- 


-2 + n +r 2 )/(- 


2 + 5ri - 3r 2 + r 2 + 3rj) , (3/2) [n + r 2 - A^/2(n + r 2 


- l)/(3e) 


M.3m 


6r 2 (- 


-2 + n +r 2 )/(- 


2 + 5ri - 3r 2 + r 2 + 3rl) , (3/2) [ri + r 2 + A^/2(ri + r 2 


- l)/(3e) 


M.4 




-3 + A 2 /(2e), 


(-9r 2 e 2 + 9rle 2 + 6neA 2 - A 4 )/ [ 3(2 - 4ri + 3r 2 - r|)e 2 - 


eX>] 



TABLE II: The corresponding eigenvalues for the first 7 critical points in Table [I] 



critical point. In Table [TT1 we present the eigenvalues for the first 7 critical points in Table HI For Point 
(M.lp), noting that its existence requires e = +1 (namely quintessence), it can be stable under condition 
(4— Ari+rf — ?" 2 )/(™ 4ri+3r 2 — 3r|) > and A > \f&e. For Point (M.lm), noting that its existence requires 
e = +1 (namely quintessence), it is unstable because the second eigenvalue is positive (nb. A is positive). 
For Point (M.2p), noting that its existence requires e = +1 (namely quintessence) and 1 < r\ — T2 < 2, 
it can be stable under condition r 2 /(2 — 5ri + 3rf + r 2 — 3r|) > and ri — r 2 < X\/2(ri — r 2 — l)/(3e). 
For Point (M.2m), noting that its existence requires e = +1 (namely quintessence) and 1 < r\ — r 2 < 2, 
it is unstable because the second eigenvalue is positive (nb. A is positive). For Point (M.3p), noting that 
its existence requires e = +1 (namely quintessence) and 1 < r\ + r 2 < 2, it can be stable under condition 
r 2 /(-2 + 5ri - 3rj + r 2 + 3r|) > and r 1 +r 2 < A^2(r x + r 2 - l)/(3e). For Point (M.3m), noting that 
its existence requires e = +1 (namely quintessence) and 1 < n + r 2 < 2, it is unstable because the second 
eigenvalue is positive (nb. A is positive). For Point (M.4), noting that its existence requires 1— A 2 / (6e) > 0, 
it can be stable under condition (-9r 2 e 2 + 9r|e 2 + 6rieA 2 - A 4 )/ [3(2 - Ar x + 3r 2 - r\)e 2 - eA 2 ] < 0. 
Finally, the eigenvalues of Points (M.5) and (M.6) are considerably involved, and hence we do not present 
them here. We find that they can exist and are stable in proper parameter-space [Tfj| . 

In summary, for the case with interaction Q = q(ap m +3(3H p m ), we find that there are two dark-energy- 
dominated attractors (M.lp) and (M.4), and four scaling attractors (M.2p), (M.3p), (M.5) and (M.6). 
These scaling attractors can help to alleviate the cosmological coincidence problem. In e.g. Q, it has been 
found that there is no scaling solution in the interacting phantom model with the familiar interaction 
Q = 3r)H p m in which rj is a constant. This fact shows that our new interaction Q = q(ap m + 3j3H p rn ) = 
iZaq ' 3<lHPm (nb- Eq. (1231) ) can bring new results to cosmology. 



IV. THE CASE OF Q = q(ap tot + 3/3Hp tot ) 



Here, we consider the case of Q — q(ap to t + 5/3Hp to t) given in Eq. (j6]). From Eq. ([8]), it is easy to find 
Ptot — '5H 2 / k 2 . Substituting into Eq. ©, we can finally obtain 

Q = ^£ (§/»-<»). (30) 

Substituting into Eq. (1171) . we find that the corresponding Q\ and Q 2 are given by 

Qi = | f |8 - as) , Q2 = - z Q/3 - as) . (31) 

Notice that q = s — 1 and s is given in Eq. (|18p . Then, substituting them into the autonomous system 
Eqs. (|14p — (fl~6)). we can find that there are 5 critical points and present the first 4 points in Table Hill All 
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Label Critical Point (x, y, z) 



T.lp +V0-3+r 4 )/e, 0, VI -r 3 



r 4 



T.lm -V(^3 + r 4 )/e, 0, VI - r 3 - r 4 



T.2p +x/(r 3 -r 4 )/e, 0, VI " ^3 + r 4 



T.2m -V(r 3 -r 4 )/e, 0, VI - r 3 + r 4 



TABLE III: The first 4 critical points for the case of Q = q{ceptot + 3/3Hptot) ■ 



the 4 points in Table Hill are scaling solutions because f2 m = z 2 > 0. The last Point (T.3) is considerably 
involved and hence we do not present it here, except to mention that it is also a scaling solution. Note 
that r3 and r 4 are given by 



2 -4a + 3/3 J 4 - 24a + 4a 2 + 20/3 - 12o/3 + 9/3 2 

4 + 6a 4 + 6a 

If 4 + 6a = 0, all the 4 critical points in Table Mil cannot exist. If 4 + 6a ^ 0, for convenience, we can 
regard r% and r 4 as the model-parameters, in place of a and /3. By reversing Eq. Q32[) . we can express a 
and j3 as functions of r3 and namely 



a = — - 



2(-l + 2r 3 + 3r 2 - 3rj) 2(-l + 2r 3 + 7r$ - 7r|) 



1 + 6r 3 + 9r 2 - 9r 2 



^ 1 + 6r 3 + 9r 2 - 9r 2 ' (33) 



3 -|- a'Hj — arj i -t- OT3 -f »r§ — arj 



Here, we briefly discuss the existence of the critical points. For Points (T.lp) and (T.lm), if e = — 1 
(namely phantom), we should have r^+r^ < by requiring x to be real. However, in this case Q m — z 2 > 1 
which is physically meaningless. Therefore, Points (T.lp) and (T.lm) can exist only for e = +1 (namely 
quintessence) and 1 > r% + r 4 > 0. Similarly, Points (T.2p) and (T.2m) can exist only for e = +1 (namely 
quintessence) and 1 > r% — r 4 > 0. Point (T.3) can exist in proper parameter-space [l6| . 
To study the stability of these critical points, by linearizing Q±, we obtain 



SQi - 4 

a; 



-/3 — as j (j)q — ~5xj — aqSs 



(34) 



where g = s — 1 and 5g = 5s, while s and 5s are given in Eq. (|23p . Substituting this <5Qi into Eqs. (|21j) 
and (1221) . the two eigenvalues of the coefficient matrix of Eqs. (j2"Tj) and ((22)) determine the stability of 
the critical point. In Table IW1 we present the eigenvalues for the first 4 critical points in Table Hill For 
Point (T.lp), noting that its existence requires e — +1 (namely quintessence) and 1 > r 3 +r 4 > 0, it can 
be stable under condition ^<^^±lJ^^r±^+^ < and | + - Xj^-^r < 0. For 

(l+3i-3-3r4)(r3+r 4 )(l-|-3r3+3r4)£ J — 2 2(r 3 +r 4 ) y 2(r 3 +r' 4 ) — 

Point (T.lm), noting that its existence requires e = +1 (namely quintessence) and 1 > r 3 + r 4 > 0, it is 
unstable because the second eigenvalue is positive (nb. A is positive). For Point (T.2p), noting that its 
existence requires e = +1 (namely quintessence) and 1 > r 3 — r 4 > 0, it can be stable under condition 
-6(r 3 +r 4 )(r3-r 4 ) 3 -i2r3(r3-r4)e 2 +i8 e 4 < and | + 9 , 3e ' . - X. fZ 3e , < 0. For Point (T.2m), noting 

(l+3r3-3r4)(r3-r 4 )(l+3r 3 +3r 4 )£ 2 — 2 2(r 3 -r 4 ) V 2(r 3 -r 4 ) — \ > ' o 

that its existence requires e — +1 (namely quintessence) and 1 > r 3 — r 4 > 0, it is unstable because the 
second eigenvalue is positive (nb. A is positive). Finally, the eigenvalues of Point (T.3) are considerably 
involved, and hence we do not present them here. We find that it can exist and is stable in proper 
parameter-space [l(| . 

So, in the case of Q — q(ap to t + 5/3Hp to t), there are 3 scaling attractors (T.lp), (T.2p) and (T.3). 
These scaling attractors can help to alleviate the cosmological coincidence problem. Of course, these 
scaling solutions are also different from the ones in the interacting quintessence or phantom model with 
the usual interaction Q = 3nHp to t in which 77 is a constant. Our new interaction Q = q(aptot + 3/377p tot ) 
brings new results. 
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Point 


Eigenvalues 
















T.lp 


-6(r 3 -r4)(r3+r4) 3 -12r3(r3+r 4 )e 2 + 186 4 


3 


+ 


3e 2 




1 3e 


(l + 3r 3 -3r4)(r3+r4)(l + 3r3+3r 4 )e a ' 


2 


2(r 3 +r 4 ) 


/ 2(r 3 +r 4 ) 


T.lm 


-6(r3-r 4 )(r3+r 4 ) 3 -12r3(r3+r 4 )e 2 + 186 4 


3 


+ 


3e 2 


+\ 


/ 36 


(l + 3r 3 -3r 4 )(r3+r 4 )(l + 3r 3 +3r 4 )6^ ' 


2 


2(r 3 +r 4 ) 


' 2(r 3 +r 4 ) 


T.2p 


-6(r3+r 4 )(r3-r4) 3 -12r3(r3-r4)e 2 +18c 4 


3 


+ 


3e 2 


-\ 


/ 3, 


(l + 3r 3 -3r 4 )(r3-r4)(l + 3r3 + 3r 4 )e i! ' 


2 


2(r 3 -r 4 ) 


/ 2(r 3 -r 4 ) 


T.2m 


-6(r3+r 4 )(r3-r 4 ) 3 -12r 3 (r3-r 4 )e 2 + 18£ 4 


3 


+ 


3e 2 


+\ 


/ 3C 


(l + 3r3-3r 4 )(r 3 -r 4 )(l + 3r 3 + 3r 4 )e a ' 


2 


2(r 3 -r 4 ) 


1 2(r 3 -r 4 ) 



TABLE IV: The corresponding eigenvalues for the first 4 critical points in Table [TTT1 



V. THE CASE OF Q = q(ap de + 3f3Hp de 



In this section, we consider the case of Q 
Eq. ([2]), one can find that 



q{apde + 3/3-ffpde) given in Eq. ((?]). Substituting it into 



-3H 



Pile 



1 + aq 

Then, substituting into Eq. ([7]), we can finally obtain 

3Hq 



(Pde + Pde + PqPde) 



Q 



1 + aq 



[PPde ~ a(p de +Pde) , 



(35) 



(36) 



which is valid for any dark energy. In the present work, the role of dark energy is played by quintessence 
or phantom. Substituting Eqs. (fTO)) and (TTTT) into Eq. (|36|) . we have 



Q = 



3Hq 



t(t> 2 + /3V 



1 + aq 

Substituting into Eq. (|1T[) . we find that the corresponding Q\ and Q2 are given by 



Qi 



3q 



2(1 + aq) 



[p — 2a) ex H 



Q-2 



3q 



2(l + aq)z 



[((3 - 2a) ex 2 + /3y 2 



(37) 



(38) 



Notice that q = s — 1 and s is given in Eq. (|18p . Then, substituting them into the autonomous system 
Eqs. (fT4"]l — (|16l) . we can find that there are 8 critical points and present the first 4 points in Table [V] All 
the 4 points in Table fVl are scaling solutions because fl m — z 1 > 0. The last 4 Points (D.3), (D.4) (D.5) 
and (D.6) are considerably involved and hence we do not present them here, except to mention that they 
are also scaling solutions. Note that r$ and re are given by 



''5 



2 + 4a -3(3 



'6 



_ v^c* 2 + 4a(10 - 3/3) + (2 - 3/3) 2 



(39) 



6a 6a 

If a = 0, all the 4 critical points in Table IVl cannot exist. If a 7^ 0, for convenience, we can regard r$ 
and r6 as the model-parameters, in place of a and /?. By reversing Eq. ([39]), we can express a and (3 as 
functions of rs and rg, namely 



9rl 



9r 2 



6r 5 - 1 



/8 = 



2(3? 



3r§ 



6r 5 + 5) 



9r§ 



9r| + 6r 5 



(40) 



Here, we briefly discuss the existence of the critical points. For Points (D.lp) and (D.lm), if e = — 1 
(namely phantom), we should have rg < r§ by requiring x to be real. However, in this case fl m = z 1 > 1 
which is physically meaningless. Therefore, Points (D.lp) and (D.lm) can exist only for e = +1 (namely 
quintessence) and rg > r$. Similarly, Points (D.2p) and (D.2m) can exist only for e = +1 (namely 
quintessence) and r§ + r§ < 0. Points (D.3), (D.4), (D.5) and (D.6) can exist in proper parameter- 
space [i~6|]. 
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Label 


Critical Point (x, 


y, ») 








D.lp 


+V(r e -r s )/e, 0, 


VI + r S 


- r 6 


D.lm 


— \/(r6-r 5 )/e, 0, 


VI + r 5 


- r 6 


D.2p 


+v/-(r 5 + r 6 )/e, 0, 


VI + r 5 


+ r e 


D.2m 


-^/-(rs + r 6 )/e, 0, 


VI + r s 


+ r 6 



TABLE V: The first 4 critical points for the case of Q — q(apde + 3/3Hpd e ) ■ 



To study the stability of these critical points, by linearizing Q±, we obtain 



SQi 



2(1 + aq) 



((3 — 2a) eSx — -^-5x H — ~~5y 

X z X 



2a)ex+^ 
x 



5q 



1 + aq 



(41) 



where q = s — 1 and 5q = 5s, while s and 5s are given in Eq. (|23p . Substituting this 5Q\ into Eqs. (|21j) 
and (1221) . the two eigenvalues of the coefficient matrix of Eqs. (j2Tj) and ((22|) determine the stability of 
the critical point. In Table EH we present the eigenvalues for the first 4 critical points in Table [Vj For 
Point (D.lp), noting that its existence requires e = +1 (namely quintessence) and r§ > rs, it can be 
stable under condition r 6 /[(3r 5 — 3r 6 — 1)(1 + r 5 + r 6 )] > and 1 — r 5 + r 6 — Xy/2(r 6 — r 5 )/(3e) < 0. 
For Point (D.lm), noting that its existence requires e = +1 (namely quintessence) and rg > rs, it is 
unstable because the second eigenvalue is positive (nb. A is positive). For Point (D.2p), noting that 
its existence requires e = +1 (namely quintessence) and r$ + r$ < 0, it can be stable under condition 
r 6 /[(l - 3r 5 - 3r 6 )(l + r 5 - r 6 )] > and 1 - r 5 - r 6 - A v /-2(r 5 + r 6 )/(3e) < 0. For Point (D.2m), noting 
that its existence requires e = +1 (namely quintessence) and r^+r^ < 0, it is unstable because the second 
eigenvalue is positive (nb. A is positive). Finally, the eigenvalues of Points (D.3), (D.4), (D.5) and (D.6) 
are considerably involved, and hence we do not present them here. We find that they can exist and are 
stable in proper parameter-space [l6j . 

So, in the case of Q = q(apd e + 3f3Hpd e ), there are 6 scaling attractors (D.lp), (D.2p), (D.3), (D.4), 
(D.5) and (D.6). These scaling attractors can help to alleviate the cosmological coincidence problem. 
Of course, these scaling solutions are also different from the ones in the interacting quintessence or 
phantom model with the usual interaction Q = Sr/Hpde in which 77 is a constant. Our new interaction 
Q = q(apd e + 3/3Hpde) brings new results. 



VI. CONCLUSION 



In the present work, motivated by the recent work of Cai and Su [14j |. we proposed a new type of 
interaction in dark sector, which can change its sign when our universe changes from deceleration to 
acceleration. We considered the cosmological evolution of quintessence and phantom with this type of 
interaction. We found that there are some scaling attractors which can help to alleviate the cosmological 
coincidence problem. Our results also showed that this new type of interaction can bring new features to 
cosmology. 
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Point 



Eigenvalues 



D.lp 
D.lm 
D.2p 
D.2m 



24r 6 (r 5 -r 6 ) 



(3r B -3r 6 -l)(l + r 5 +re) ' 

24r 6 (r 5 -r 6 ) 

(3r 5 -3r 6 -l)(l + r 5 +r 6 ) ' 

2ir 6 (r 5 +r 6 ) 

Cl-3r 6 -3r 6 )(l+r 5 -r 6 ) ' 

24r 6 (r 5 +r 6 ) 

(l-3r 5 -3r 6 )(l + r 5 -r e ) ' 



I - r 5 + r 6 - A^/2(r 6 - r 5 )/(3e) 
I - r 5 + r 6 + A ^2(7-6 - r 5 )/(3e) 



(3/2) [ 
(3/2) [ 

(3/2) [l - r s - r 6 - A N /-2(r 5 + r 6 )/(3e 
(3/2) [l - r 5 - r 6 + A^-2(r 5 + r 6 )/(3e 
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